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We have built an effective Walecka-type hadronic Lagrangian in which the hadron 
^ ' masses and the density dependence of the coupling constants are deduced from the 

quark dynamics using a Nambu-Jona-Lasinio model. In order to stabilize nuclear 
matter an eight-quark term has been included. The parameters of this Nambu-Jona- 
Lasinio model have been determined using the meson properties in the vacuum but 
also in the medium through the omega meson mass in nuclei measured by the TAPS 



I collaboration. Realistic properties of nuclear matter have been obtained. 
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1 Introduction 

One of the most fascinating challenges of nuclear physics is the description of 
nuclear matter and nuclei starting from Quantum ChromoDynamics (QCD). 
Even if important progress have been made in QCD calculations on the lat- 
tice, such a description is not yet available. Models incorporating the most 
prominent features of QCD have to be used. 

Quark models, like, for example, the quark-meson coupling model[l,2] which 
describes nuclear matter as non overlapping nucleon MIT bags interacting 
through the self-consistent exchange of meson fields directly coupled to the 
quarks, or the model of Bentz and Thomas [3] where the nucleons are gener- 
ated as quark-diquark states in a Nambu-Jona-Lasinio (NJL) model with an 
infrared cutoff to simulate the confinement, are gaining more and more success 
in shedding some light on the dynamical behaviour of nucleons in nuclei. 
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On the other hand, starting directly at the hadronic level, Kohn-Sham-like 
density functional theories with constraints on the energy functional coming 
from low energy QCD, like that developed by Finelli, Kaiser, Vretenar and 
Weise[4], lead to an accurate description of nuclear matter and finite nuclei. 

Another possibility is to apply the strategy of effective field theories where 
low energy effective hadronic Lagrangians are obtained by integrating out the 
degrees of freedom lying above the energy scale considered. This decimation 
leads to density-dependent masses and couplings in the hadronic Lagrangian 
[5]. This density-dependence should reflect essentially the underlying dynam- 
ics of the quarks. Such a Lagrangian with density dependent masses and cou- 
pling constants determined according to Brown and Rho scaling[6] has been 
proposed by Brown, Song, Min and Rho [7]. The calculation reported in [7], 
assuming a scaling law leading to a decreasing of the vector meson mass of 
approximately 20% at saturation, enables a realistic description of bulk prop- 
erties of nuclear matter. 

Recently, new experimental results for the in-medium uj meson mass would 
suggest a small decreasing of approximately 10-15% at saturation density. In 
particular, the TAPS collaboration [8] found an in-medium mass of m* = 
722i4(stat)^g 5 (syst) MeV at 0.6 times the saturation density (which is com- 
patible with a 14% dropping at saturation when a linear density-dependence 
is assumed) in photoproduction experiments and Naruki et al. [9] found a 9% 
decrease of the in-medium u mass at saturation in 12 GeV proton-nucleus 
reactions. 

In this paper, we explore the possibility of obtaining a realistic description 
of bulk properties of nuclear matter in the model of Brown, Song, Min and 
Rho [7] with a density dependence of the in-medium u mass in accordance 
with recent experimental indications. In addition, in order to allow deviations 
from the Brown and Rho scaling, the hadron masses and the density depen- 
dence of the coupling constants have been deduced directly from the quark 
dynamics using an in-medium NJL model. We have chosen here to use the 
NJL model [10] since, despite some shortcomings like the lack of confinement, 
with its dynamical quark mass generation and the in-medium chiral symme- 
try restoration, it allows to take into account an important part of quark 
dynamics. 

The paper is organised as follows. In section 2, we recall the formalism of the 
effective hadronic model with density-dependent couplings and masses used 
to determine bulk properties of infinite nuclear matter. The section 3 will 
be devoted to the determination of density dependences for hadron masses 
and meson couplings from a NJL model. We present and discuss numerical 
results in section 4, where some nuclear matter bulk properties (saturation 
point, incompressibility, nucleon effective mass ...) are presented. We conclude 
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in section 5. 



2 Effective hadronic mean-field model with density-dependent pa- 
rameters 

We use an effective hadronic Lagrangian similar to that used by Brown, Song, 
Min and Rho [7]. It is formally identical to the Walecka onefll] but with 
density-dependent masses and couplings, which reflect the underlying physics 
induced by the quark dynamics: 

C had = ^j[-f^{id^ - gujNN^V^) - (M* N - g aNN a* a <f))} ip 

+\{d^<f> - mfcf) + \rrtfV,V» - \f^\ (1) 

where ip, <fi and V 1 represent respectively the nucleon, a and u> meson fields 
and, as usual, F^ v = d^V u — d u V^. This Lagrangian contains only two free- 
parameters, QaNN an d g W NN, which will be adjusted at the end of the calcula- 
tion to reproduce the position of the saturation point. 

As already indicated, the density-dependent a* functions and mass parame- 
ters, Mjy-, m*, m*, would contain information coming from the quark level. 
The problem is now to determine these quantities in such a manner they re- 
flect the underlying quark dynamics as better as possible. If a fully satisfactory 
description of nucleons and nuclear matter from quark and gluon fields would 
be available, this stage of calculation would be unnecessary since the fields 
and parameters entering Chad could be expressed directly from the quark level 
ones. However, it is not actually the case and approximations have still to be 
made. For example, in their work, Brown, Song, Min and Rho [7] have pre- 
supposed a density dependence of the Lagrangian parameters in accordance 
with Brown and Rho scaling. In the present exploratory work, we have cho- 
sen to assume for the values of the masses and couplings the same baryonic 
density dependence as that which would be obtained in quark matter for the 
same density. As explained in the next section, we have used for this job a 
NJL model. Obviously, in the real world, quark matter cannot exist at low 
temperatures and densities but we expect that, despite some shortcomings, 
the approximation used is not so bad since we have verified that, for example, 
the density dependence of the quark condensate obtained is very similar to 
that obtained in ref[3] using a quark-diquark description of the nucleons in 
nuclear matter. 

In the rest frame of symmetric nuclear matter, the equation for the nucleon 
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field in mean-field approximation is: 



- g U NNOtil V o + 7 s o - (M* N - g aNN al<po) i/j = 0, 



(2) 



where the scalar and vector classical fields 0o and Vq can be obtained from 
the Lagrange equations : 



</>0 
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(3) 
(4) 



with the scalar and baryonic nuclear densities, p s and ps, defined as usual as 

Ps = (^) , (5) 
Pb = UH) ■ (6) 



For each value of the baryonic density, p s will be determined self-consistently 
by minimizing the energy density £ (Eq.10) with respect to p s . In Eq.2, S 
represents the rearrangement term given by : 

v _ 91nn<P 2 b d(K/m*J 8M* N g 2 aNN afp 2 s d(m*Ja* a ) 

^0 — — Ps— =^3 — > V<) 



m: 



dpi 



dp 
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expression in which a*, a*, M^, m* and m* are functions of ps- As in the 
standard Walecka model[ll], Eq.2 is equivalent to the Dirac equation for a 
nucleon of mass : 



m* N = M* N - g aN NaU 



(8) 



Note that the vector potential Uq is given by : 
Uq = gu;NN&tyo — So- 



(9) 



As usual, the energy-density £, pressure p and incompressibility parameter 
K can be obtained from the energy-momentum tensor, providing, for Fermi 
momentum hp '■ 



c_ 91nn< 2 P 2 b , 91nnK 2 P 2 s , 7 
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p = 



91nn< 2 P 2 b _ 91nn< 2 P 2 s 
2m* 2 2m* 2 



- Pb^o 



(11) 
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kpe* F (2kp — 3m^) + 3m^ In 



( 




) 



d 2 S 



(12) 



where e* F = ykp + m*^ and 7 = 4 in symmetric nuclear matter. As already 
proved by Brown, Song, Min and Rho [7], the thermodynamical consistency 
holds, whatever the density dependence used in Eq.l, as long as the rearrange- 
ment terms (Eq.7) are not forgotten. 



3 Quark dynamics in a Nambu Jona-Lasinio model 

3.1 In medium quark mass 

We consider the following chirally invariant two-flavor NJL Lagrangian[10] : 



where q denotes the quark field with two flavor (Nf = 2) and three color 
(N c = 3) degrees of freedom and m is the diagonal matrix of the current 
quark masses (here in the isospin symmetric case). The second and third 
terms of Eq.13 represent local four-quark interactions while the later one is an 
eight-quark interaction. As explained in Section 4, the eight-quark term, while 
contributing very weakly to the energy, is essential here because its influence on 
the density dependences of the hadronic Lagrangian (Eq.l) enables to stabilize 
nuclear matter. In fact, this term provides an additional density dependence of 
the Lagrangian parameters which is important for the scalar polarizability of 
the nucleon. It is interesting to note that such a stabilizing effect has already 
been reported by Bentz and Thomas [3] in a quite different context, and by 
Mishustin et al. in a nucleonic NJL model [12]. Other high-order interaction 
terms could be added to the Lagrangian 13. The effect of some of them is 
actually under consideration but, in this exploratory work, only that appearing 
in 13 will be taken into account. 



£n.jl = q [il^ -m ]q + g x [(qq) 2 + (gry 5 rg) 2 J - g 2 (qi^q) 2 
+93 \(qq) 2 + (qil5rq) 2 ] (qi^q) 2 , 



(13) 
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The Dirac equation for a quark in mean-field approximation is given by : 

i^d» -mo- 2(? 2 7 (g 7 °g) + 2 9l (qq) + 2g 3 (qq) (g 7o g) 2 l q = 0, (14) 



which defines a dynamical constituent-quark mass : 

m = m - 25(1 11 H ^— I (gg) = m - g 1 {qq} , (15) 



generated by a strong scalar interaction of the quark with the QCD vacuum. 
In the gap equation (Eq.15), the quark condensate (qq) can be written as : 

(qq) =-iJ ^JrS(k), (16) 



where here Tr denotes traces over color, flavor and spin. In Eq.16, S(k) rep- 
resents the in-medium quark propagator defined as : 

S(k) = — + i7r ^ k Z + m S (ho - E k ) 9 (k F - Ikl) , (17) 



where k* = k„ - 2(g 2 - g 3 (qq)) (q lfl q), E* k = Vk* 2 + m 2 , E k = E£ + 2(g 2 - 
93 (QQ)) (QloQ) an d kp is the quark Fermi momentum. It has to be noted that we 
are working in quark matter in the NJL model. The baryonic density is related 
to the total quark density by p# = The quark condensate is divergent due 
to the loop integrals and requires an appropriate regularization procedure. As 
many authors[13,14], we introduce a three-momentum cutoff A which has the 
least impact on medium parts of the regularized integrals, in particular at zero 
temperature [14]. In fact, since the model is non renormalisable, the cut-off A 
is just an additional parameter. 

Thus, after the regularization procedure, the quark condensate is given at each 
density by : 

, . N f N c A f mk 2 dk , , 

<^ = -^/— < 18 > 

kp 



self-consistent equation which is equivalent to a minimization of the energy 
density of quark matter with respect to the quark condensate. 



6 



3.2 Mesons masses and couplings 



As usual, the mesons will be obtained by solving the Bethe-Salpeter equation 
in the quark- ant iquark channel. Since this is the standard procedure, we will 
only list the results which are needed later on. First, we define the quark- 
antiquark polarization operator in the M = n,a,u channel by : 

n M (g 2 ) = -ij -^fTr [T M iS{p + q/2)T M iS{p - q/2)} , (19) 



where the vertex Tm = il5i~ k , il and 27^ stand for respectively the pion, sigma 
and omega mesons. Note that, in the vector channel the Lorentz structure of 
the polarization operator is H% v (q 2 ) = [—g^ + ^f-j Yi^q 2 ). The in-medium 
meson masses, m* M , and meson-quark-quark coupling constants, g Mqq , are then 
determined by the pole structure of the T-matrix, i.e. by the conditions: 



1 - K M U M (q 2 = m M ) = 0, 



(20) 



9*Mqq 



<ffl M (g 2 

dq 2 



(21) 



a 2 =m *2 
1 " l M 



where Km = gi, <?i and 2g 2 — 2g 3 (qq) 2 respectively for the ir, a and uj mesons. 
Note that for the determination of the polarization operator (Eq.19), we have 
used the same regularization procedure as for the quark condensate. 



In the scalar channel, the polarization is given by: 



n CT (g 2 ) = 



(qq) 



m 



+ N c N f (q 2 -4m 2 )I(q 2 ), 



(22) 



with 

- & T '^f-^^- <23) 

4(A 2 +m 2 ) 



In the vector channel, one obtains: 



4(pf,+m 2 ) 



IWH-A^ / ^(l + ^l-^h.. (24) 

4(A 2 +m 2 ) 
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We also need the pion mass and decay constant to adjust the model parame- 
ters. In the pseudo-scalar channel, the polarization reads: 



U 7r (q 2 ) = ^ + N c N f q 2 I(q 2 ), (25) 



and the pion decay constant is : 

U = N c N f g nqq mI(q 2 = ml). (26) 

3. 3 Density dependences for the hadronic model 



We now have to define how the density dependences of the hadronic La- 
grangian parameters (Eq. 1) are deduced from quantities calculated in NJL 
model. First, for every value of the baryonic density, the u> and a meson 
masses m* , m* are taken directly from the NJL calculation (Eq. 20) in which 
the value of the quark condensate has been obtained solving Eq.18. Second, 
we have to determine M^. Even if progress have been made (see for example 
[3,15]), a realistic description of the nucleon in a NJL model is still an open 
problem. Here, we have assumed that is directly related to the quark 
condensate, with the same relation as found in finite-density QCD sum-rule 
calculations [16]. To leading order, which should be valid at densities below 
and around the saturation density of nuclear matter, one has then : 

Mn = M_ (07) 
M N (qq) ' 1 ] 



where (qq) is the value of the quark condensate at baryonic density ps obtained 
by solving Eq.18, (qq) is the same quantity in vacuum and M N = 939 MeV 
is the free nucleon mass. 

It was predicted in [5] that the dependence of the vector-meson mass on the 
quark condensate should change around saturation from approximately (qq) 5 
below saturation to (qq) above. The dynamical processes leading to this change 
would probably also affect the nucleon mass at the same density. However, 
since in this work we are interested only in densities below and around satura- 
tion, we have assumed that Eq. 27 remains valid at every density considered. 
Third, in a simple N c counting in a naive quark model, the quark-meson and 
nucleon-meson couplings in vacuum are proportional, as argued in [13]. At low 
densities, it seems quite reasonable to assume that the same relation holds 
approximately. Such a proportionality holds also, for example, in the QMC 
model[l,2]. The a* and a* functions entering in the hadronic Lagrangian 



8 



would thus be taken as: 




Qcrqq 



(28) 




Qwqq 



(29) 



where g* , g^ qq and their values in vacuum g aqq , g^qq, are obtained from Eq.21 
with the value of the quark condensate which minimizes the energy density of 
quark matter. 



4 Results and discussion 

4-1 NJL model parameters 

At the quark level, we have to fix five parameters for the NJL model : the cut- 
off A, the bare quark mass m , and the coupling constants gi, g 2 and g 3 . For a 
given value of the cutoff A, we need four constraints. As usual, the pion mass 
m n = 135 MeV, the pion decay constant f w = 92.4 MeV and the uo meson mass 
= 782 MeV in vacuum are used. Since the eight-quark term has been in- 
troduced in order to stabilize nuclear matter, at least one in-medium physical 
quantity should be constrained in the fitting procedure. We have chosen here 
to take into account the recent result obtained by the TAPS collaboration [8] 
for the ijj meson mass in nuclei, m* (p B = 0.6po) — 722l 4 (stat)l5 5 (syst) MeV 
(where po is the saturation density of nuclear matter). Thus, the parameters 
m o, 9i, 92 and g 3 are obtained by fitting the pion and the uj meson properties 
in vacuum and the in-medium uj meson mass m* (p B = 0.6p ) = 722 MeV. 
We have verified that the results are qualitatively and quantitatively similar 
if we choose another constraint on the in-medium u mass in accordance with 
the recent experimental data of Naruki et al.[9] on 12 GeV proton- nucleus 
reactions. With these constraints, the cutoff A cannot be unambiguously de- 
termined and, as often done, we have considered here several different values. 



4-2 Hadronic constraints and observables 

On the other hand, at the hadronic level, g a NN and g^NN have been determined 
to reproduce the empirical saturation point of nuclear matter, i.e., £ / Pb — 
M N = -16 MeV and p = at p B = 0.17 fm" 3 . 
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For each value of A considered, in order to probe the description of nuclear 
matter obtained, we have calculated several physical quantities: the effective 
nucleon mass m* N , the incompressibility parameter K and the slope of the real 
part of the energy dependence of the nucleon-nucleus optical potential U /M N 
at saturation. For the compressibility of nuclear matter, the main information 
comes from the energies of the monopolar modes over many nuclei. The com- 
monly accepted value K = 210 ± 30 MeV, has been deduced by Blaizot[17], 
in the framework of non-relativistic Hartree-Fock RPA models with effective 
interactions. In a relativistic framework, it has been shown[18] that some ef- 
fective Lagrangians leading to higher values of the incompressibility (close 
to 300 MeV) could predict correctly the isoscalar giant resonance energies in 
medium and heavy nuclei. Thus, we can estimate that a realistic value of the 
incompressibility parameter should be in the interval: K = 250 ± 50 MeV. 
Concerning the effective nucleon mass, Furnstahl, Rusnak and Serot[19] have 
shown that for models without isoscalar tensor coupling, there is a strong cor- 
relation between the value of the nucleon effective mass at saturation density 
and the spin-orbit splitting in nuclei. They found that, to accurately repro- 
duce the empirical splittings, one requires a value of m* N /M N between 0.58 
and 0.64 at saturation density. Finally, the slope of the energy dependence of 
the real part of the nucleon-nucleus optical potential, which in a relativistic 
mean-field approximation is Uq/Mn, has been extracted [20, 2 1,22] from exper- 
imental data[23] and has been found to be equal to 0.30 in Ref. [20,21] and 0.35 
in Ref. [22], up to 100 MeV of incident kinetic energy, at saturation. Thus, it 
is reasonable to expect that realistic values of U /M N should lie between 0.25 
and 0.40. 



4-3 Numerical results 

Each value of A corresponds to a value of the constituent quark mass in vac- 
uum, and we have thus chosen to display the results as functions of m instead 
of A. With the constraints chosen here at the quark level, the saturation point 
cannot be reproduced for m < 460 MeV. We have plotted on Fig.l the dimen- 
sionless nucleon effective mass m* N /M N , the incompressibility parameter K 
and the slope of the energy dependence of the real part of the nucleon-nucleus 
optical potential U /M N at saturation as functions of m. 

The shaded areas correspond to the bounds on the empirical values. As we can 
see, for m « 465 - 470 MeV (A = 572 ± 1 MeV), the three physical quantities 
m* N /M N) K and U /M N are all in good agreement with the empirical values, 
and keep reasonable values at least up to m 500 MeV. 

As an example, we have chosen to consider the central value m = 467 MeV 
represented by the vertical dashed line in Fig.l. The parameters of the NJL 
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Fig. 1. Dimensionless effective nucleon mass m^/M^, incompressibility parameter 
K and slope of the energy dependence of the nucleon-nucleus optical potential 
Uq/Mn shown as functions of the constituent quark mass m. The shaded areas 
correspond to the bounds on the empirical values (see text) and the vertical dashed 
line indicates the central value m = 467 MeV. 

model are then A = 572 MeV, m = 5.6 MeV, g l = 8.22 GeV" 2 , g 2 = 21.4 
GeV" 2 and g 3 = 1.08 10 4 GeV" 8 . At the hadronic level, the values of the 
coupling constants are g a NN = 6.3 and g U NN — 13. The saturation curve for 
symmetric nuclear matter is shown on Fig.2. 

The values obtained for the effective nucleon mass, the incompressibility pa- 
rameter and the slope of the optical potential at saturation are m* N /MN = 0.6, 
K = 221 MeV, and U /M N = 0.32. Some remarks are in order. 

First, the value of the g 3 parameter is such that g 3 N^ p\ j '4g 1 = 0.02 pro- 
viding an eight-quark term which contributes only for 2% of the constituent 
quark mass at saturation (see Eq.15). However, this contribution is essential 
for stabilizing the system. With g 3 = 0, nuclear matter would saturate with 
an unrealistically weak effective nucleon mass at saturation density (~ 30% of 
the free nucleon mass). The eight quark term tends to smooth the curvature 
of the quark condensate (and that of in this model) with respect to den- 
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Fig. 2. Saturation curve of nuclear matter obtained for m = 467 MeV. 

sity. The sensitivity of the model to this term can be understood as follows. 
The saturation is resulting as a compensation between the scalar, vector and 
rearrangement term (Eq. 7) contributions in pressure, and the rearrangement 
term is very sensitive to the derivatives of the masses and couplings with re- 
spect to density. The additional density dependence produces a variation of 
the Eo term of about 40% at saturation. So, the difference between g% = 
and #3 lies principally in the modification of the derivatives of the masses 
and couplings rather than of these functions themselves. 

Second, at the quark level, the model provides a quark condensate in the 

1 /3 

vacuum (uu) = —240 MeV in good agreement with the lattice calculations: 
(uu) 1/3 = -(231 ± 4 ± 8 ± 6) MeV[24]. 

Third, the constituent quark mass in vacuum is relatively large. This value 
is appreciably larger than one third of the free nucleon mass, but since this 
quantity is not an observable, it should be model-dependent. Note that such 
large values are commonly considered in recent NJL calculations (see for ex- 
ample [14,25]). In particular, a large value has also been obtained[25] in a 
NJL model including meson loop corrections generated via a 1/N C expansion 
of the self-energy in the next to leading order, with parameters constrained by 
the experimental data on the electromagnetic form factor F n in the time-like 
region. Such a large mass doesn't prevent from obtaining a good description 
of the baryonic spectrum, for example using a quark-diquark model as done 



12 



in ref [26,27]. Moreover, a large quark mass prevents the omega meson to be 
unstable against decay into a quark-antiquark pair since m* is always lower 
than 2m for each density. 

Fourth, one can see on Fig. 2 that the saturation curve is rather hard at den- 
sities above po- This can be due to the fact that M^, taken here proportional 
to the quark condensate, is decreasing too fast at high densities. This could be 
compensated, for example, by introducing perturbatively higher order terms 
in the scalar and vector fields, as done in [7] . 



5 Conclusion 

We have investigated the properties of nuclear matter in a relativistic mean 
field model with density-dependent masses and couplings, similar to that used 
in [7], in which we have replaced the Brown and Rho scaling by a direct 
calculation of meson masses and couplings in a NJL quark model. The NJL 
model including four and eight quark interaction terms has been constrained 
to reproduce the recent TAPS result for the in medium uj meson mass in 
addition to the vacuum pion and u meson properties. At the hadronic level, 
the two free parameters have been fixed to reproduce the empirical saturation 
point. 

At the quark level, the eight quark term is essential for obtaining realistic 
saturation properties. At saturation, the nucleon effective mass, the incom- 
pressibility parameter and the slope of the energy dependence of the nucleon- 
nucleus optical potential obtained are all in good agreement with the empirical 
data. 

Even if more work is needed in this direction and in a more fundamental one, 
this result is a very encouraging one since with only a few free-parameters, a re- 
alistic description of saturation properties of nuclear matter has been obtained 
from a hadronic Lagrangian constrained by a quark model which reproduces 
vacuum and in-medium meson properties. 
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